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Abstract 

Several approaches to testing the hypothesis that two histograms are drawn 
from the same distribution are investigated. We note that single-sample continuous 
distribution tests may be adapted to this two-sample grouped data situation. The 
difficulty of not having a fully-specified null hypothesis is an important consideration 
in the general case, and care is required in estimating probabilities with "toy" Monte 
Carlo simulations. The performance of several common tests is compared; no single 
test performs best in all situations. 

1. Introduction 

Sometimes we have two histograms and are faced with the question: "Are they 
consistent?" That is, are our two histograms consistent with having been sampled 
from the same parent distribution. For example, we might have a kinematic distri- 
bution in two similar channels that we think should be consistent, and wish to test 
this hypothesis. Each histogram represents a sampling from a multivariate Poisson 
distribution. The question is whether the means are bin-by-bin equal between the 
two distributions. Or, if we are only interested in "shape", are the means related 
by the same scale factor for all bins? We investigate this question in the context of 
frequency statistics. 

For example, consider Fig. 1. Are the two histograms consistent or can we 
conclude that they are drawn from different distributions? 

* This work supported in part by the U.S. Department of Energy under grant DE-FG02- 
92-ER40701. 
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Fig. 1. Two histograms (blue squares and red circles) to be compared for consistency. 



There are at least two variants of interest to this question: 

1. We wish to test the hypothesis: 

Hq: The means of the two histograms are bin-by-bin equal, against 
Hi: The means of the two histograms are not bin-by-bin equal. 

2. We wish to test the hypothesis: 

Hq. The densities of the two histograms are bin-by-bin equal, against 

H[: The densities of the two histograms are not bin-by-bin equal. 

In the second case, the relative normalization of the two histograms is not an issue: 
we only compare the shapes. 

It may be noted that there are a large variety of tests that attempt to answer the 
question of whether a given dataset is consistent with having been drawn from some 
specified continuous distribution. These tests may typically be adapted to address 
the question of whether two datasets have been drawn from the same continuous 
distribution, often referred to as "two-sample" tests. These tests may further be 
adapted to the present problem, that of determining whether two histograms have 
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the same shape. This situation is also referred to as comparing whether two (or 
more) rows of a "table" are consistent. The datasets of this form are also referred to 
as "grouped data". 

Although we keep the discussion focussed on the comparison of two histograms, 
it is worth remarking that many of the observations apply also to other situations, 
such as the comparison of a histogram with a model prediction. 



2. Notation 

Wc assume that we have formed our two histograms with the same number of 
bins, k, with identical bin boundaries. The bin contents of the "first" histogram are 
given by realization u of random variable U, and of the second by realization v of 
random variable V. Thus, the sampling distributions are: 

1=1 

k Vi 



where the vectors and v are the mean bin contents of the respective histograms. 
We define: 

k 



Nu = Ui, total contents of first histogram, 
1=1 

k 

Nv = Vi, total contents of second histogram, 

i=l 

k 



i=l 
k 

VT = {Ny) ^^Pi, 
1=1 

ti = Ui -\- Vi, i — \, . . . , k. 



We are interested in the power of a test, at any given confidence level. The power 
is the probabihty that the null hypothesis is rejected when it is false. Of course, the 



power depends on the true sampling distribution. In other words, the power is one 
minus the probabihty of a Type II error. The confidence level is the probability 
that the null hypothesis is accepted, if the null hypothesis is correct. Thus, the 
confidence level is one minus the probability of a Type 1 error. In physics, we usually 
don't specify the confidence level of a test in advance, at least not formally. Instead, 
we quote the P-value for our result. This is the probability, under the null hypothesis, 
of obtaining a result as "bad" or worse than our observed value. This would be the 
probability of a Type I error if our observation were used to define the critical region 
of the test. 

Note that wc arc dealing with discrete distributions here, and exact statements of 
frequency arc problematic, though not impossible. Instead of attempting to construct 
exact statements, our treatment of the discreteness will be such as to err on the 
"conservative" side. By "conservative" , we mean that we will tend to accept the null 
hypothesis with greater than the stated probabitlity. It is important to understand 
that this is not always the "conservative" direction, for example it could mislead us 
into accepting a model when it should be rejected. 

We will drop the distinction between the random variable (upper case symbols 
U and V) and a realization (lower case u and v) in the following, but will point out 
where this informality may yield confusion. 

The computations in this note are carried out in the framework of the R statistics 
package [1]. 

2.1 Large Statistics Case 

If all of the bin contents of both histograms are large, then we may use the 
approximation that the bin contents are normally distributed. 

Under Hq, 

(ui) = (vi) = Hi, i = 1, . . . , A;. 

More properly, it is (Ui) — /li, etc., but we are permitting Ui to stand for the random 
variable as well as its realization, as noted above. Let the difference in the contents 
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of bin i between the two histograms be: 



Ai = Ui- 



and let the standard deviation for Aj be denoted (Tj. Then the samphng distribution 

of the difference between the two histograms is: 



P(A) = 




This suggests the test statistic: 



If the (Tj were known, this would simply be distributed according to the chi-square 
distribution with k degrees of freedom. The maximum-likelihood estimator for the 
mean of a Poisson is just the sampled number. The mean of the Poisson is also its 
variance, and we will use the sampled number also as the estimate of the variance in 
the normal approximation. 

We suggest the following algorithm for this test: 

1. For af form the estimate 



li Ui — Vi — for bin i, the contribution to the sum from that bin is zero. 

3. Estimate the P-value according to a chi-square with k degrees of freedom. Note 
that this is not an exact result. 

If it is desired to only compare shapes, then the suggested algorithm is to scale 
both histogram bin contents: 





(Ui + Vi). 
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Statistic T is thus evaluated according to: 




k 
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1. Let 



N = 0.5{N^ + N^). 



Scale u and v according to: 

Ui^u'i^ Ui{N/Nu) 
Vi^v[^Vi{N/Ny). 



2. Estimate a'f with: 



..2 / f 



3. Statistic T is thus evaluated according to: 



T = y^ 



k I Jii. — JLL 



1=1 m^N^ 

3. Estimate the P-value according to a chi-square with k — 1 degrees of freedom. 
Note that this is not an exact result. 

Due to the presence of bins with small bin counts, we might not expect this 
method to be especially good for the data in Fig. 1, but we can try it anyway. 
Table I gives the results of applying this test, both including the normalization and 
only comparing shapes. 

Table 1. Results of tests for consistency of the two datasets in Fig. 1. The tests below 
the lines are described in Section 3. 



Type of test 


T 


NDOF 


P(x2 > T) 


P-value 


Absolute comparison 


29.8 


40 


0.88 


0.86 


X^ Shape comparison 


24.9 


39 


0.96 


0.95 


Likelihood Ratio Shape comparison 


25.3 


39 


0.96 


0.96 


Kolmogorov-Smirnov Shape comparison 


0.043 


39 


NA 


0.61 


Bhattacharyya Shape comparison 


0.986 


39 


NA 


0.97 


Cramer- Von- Miscs Shape comparison 


0.132 


39 


NA 


0.45 


Anderson-Darling Shape comparison 


0.849 


39 


NA 


0.45 


Likelihood value shape comparison 


79 


39 


NA 


0.91 
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In the column labeled "P- value" an attempt is made to compute (by simulation) 
a more reliable estimate of the probability, under the null hypothesis, that a value 
for T will be as large as that observed. This may be compared with the P{x^ > T) 
column, which is the probability assuming T follows a distribution with NDOF 
degress of freedom. 

Note that the absolute comparison yields slightly poorer agreement between the 
histograms than the shape comparison. The total number of counts in one dataset 
is 492; in the other it is 424. Treating these as samplings from a normal distribution 
with variances 492 and 424, we find a difference of 2.2 standard deviations or a two- 
tailed P-value of 0.025. This low probability is diluted by the bin-by-bin test. Using a 
bin-by-bin test to check whether the totals are consistent is not a powerful approach. 
In fact, the two histograms were generated with a 10% difference in expected counts. 

The evaluation by simulation of the probability under the null hypothesis is 
in fact problematic, since the null hypothesis actually isn't completely specified. 
The problem is the dependence of Poisson probabilities on the absolute numbers of 
counts. Probabilities for differences in Poisson counts are not invariant under the 

total number of counts. Unfortunately, we don't know the true mean numbers of 
counts in each bin. Thus, we must estimate these means. The procedure adopted 
here is to use the maximum likelihood estimators (see below) for the mean numbers, 
in the null hypothesis. We'll have further discussion of this procedure below - it does 
not always yield valid results. 

In our example, the probabilities estimated according to our simulation and the 
probabilities according to a distribution are close to each other. This suggests 
the possibility of using the probabilities - if we can do this, the problem that we 
haven't completely specified the null hypothesis is avoided. We offer the following 
conjecture: 

Conjecture: Let T be the test statistic described above, for either the absolute 
or the shape comparison, as desired. Let Tc be a possible value of T (perhaps the 
critical value to be used in a hypothesis test). Then, for large values of Tq. 

P{T <T,)>P{T< n\x\T, ndof)) , 

where P (T < Tc\x^{T, ndof}) is the probability that T < Tc according to a dis- 
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tribution with ndof degrees of freedom (either k or k — 1, according to which test is 
being performed). 

We'll only suggest an approach to a proof, which could presumably also be used 
to develop a formal condition for to be "large" . The conjecture also appears to 
be true anecdotally, and for interesting values of Tc, noting that it is large values of 
Tc that we care most about for an interesting hypothesis test. 

We provide some intuition for the conjecture by considering the case of one bin. 
For simplicity we'll also suppose that u = fi and that fi is small 1, say). Since 
wc arc interested in large values of the statistic, wc arc interested in the situation 
where one ofu,v is large, and the other small (since n is small). Suppose it is u that 
is large. Then 

T = ^ fa u. 

u + V 

For given v (0, say), the probability of T is thus 

This may be compared with the chi-square probability distribution for one degree of 
freedom: 

The ratio is, dropping the constants: 

P(T) /x^e^/2VT _ exp[r(i + ln/.)] 
P(T = x2) T! ~ VTT{T) 

which approaches zero for large T, for any given /i. We conclude that the conjecture 
is valid in the case of one bin, and strongly suspect that the argument generalizes to 
multiple bins. 

According to the conjecture, if we use the probabilities from a distribution 
in our test, the error that we make is in the "conservative" direction (as long as 
Tc is large). That is, we'll reject the null hypothesis less often than we would with 
the correct probability. It should be emphasized that this conjecture is independent 
of the statistics of the sample, bins with zero counts are fine. In the limit of large 
statistics, the inequality approaches equality. 



Lest we conclude that it is acceptable to just use this great simplification in all 
situations, we hasten to point out that it isn't as nice as it sounds. The problem 
is that, in low statistics situations, the power of the test according to this approach 
can be dismal. That is, we might not reject the null hypothesis in situations where 
it is obviously implausible. 

We may illustrate these considerations with some simple examples, see Fig. 2. 
The plot for high statistics on the left shows excellent agreement between the actual 
distribution and the distribution. The lower statistics plots in the middle and 
right, for two different models, show that the chi-square approximation is very con- 
servative in general. Thus, using the chi-square probability lacks power in this case, 
and is not a recommended approximation. 




60 80 100 120 140 60 80 100 120 140 60 80 100 120 140 

Tc Tc Tc 



Fig. 2. Comparison of the actual (cumulative) probability distribution for T with the 
chi-square distribution. The solid blue curves show the actual distributions, and the 
dashed red curves the chi-square distributions. All plots are for 100 bin histograms, 
(a) Each bin has mean 100. (b) Each bin has mean 1. (c) Bin j has mean 30/j. 



3. General Case 

If the bin contents are not necessarily large, then the normal approximation may 
not be good enough. There are various approaches we could take in this case. We'll 
discuss and compare several possibilities. 
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3.1 Combining Bins 

A simple approach is to combine bins until the normal approximation is good 
enough. In many cases this doesn't lose too much statistical power. It may be 
necessary to check with simulations that probability statements are valid. Figure 3 
shows the results of this approach on the data in Figure 1, as a function of the 
minimum number of events per bin. The comparison being made is for the shapes. 
The algorithm is to combine corresponding bins in both histograms until both have 
at least "minBin" counts in each bin. 



in 




minBin minBin 



Fig. 3. Left: The blue dots show the value of the test statistic T, and the red pluses 
shows the number of histogram bins for the data in Fig. 1, as a function of the 
minimum number of counts per histogram bin. Right: The P- value for consistency 
of the two datasets in Fig. 1. The red pluses show the probability for a chi-square 
distribution, and the blue circles show the probability for the actual distribution, 
with an estimated null hypothesis. 



3.2 Testing for Equal Normalization 

An alternative is to work with the Poisson distributions. Let us separate the 
problem of the shape from the problem of the overall normalization. In the case of 
testing equality of overall normization, there is a well-motivated choice for the test 
statistic, even for low statistics. 
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To test the normalization, we simply compare totals over all bins between the 
two histograms. Our distribution is 



The null hypothesis is Hq : /it — J^T, to be tested against alternative Hi : ^ ^T- 
We are thus interested in the difference between the two means; the sum is effectively 
a nuisance parameter. That is, we are interested in 




This probability now permits us to construct a uniformly most powerful test of 
our hypothesis (Ref. 2). Note that it is simply a binomial distribution, for given 
N. The uniformly most powerful property holds independently of N, although the 
probabilities cannot be computed without A^. 

The null hypothesis corresponds to jj,T = i^T, that is: 

P{Ny\Nu + Ny^N) = 




For our example, with = 916 and Ny = 424, the P-value is 0.027, assuming a 
two-tailed probability is desired. This may be compared with our earlier estimate 
of 0.025 in the normal approximation. Note that for our binomial calculation we 
have "conservatively" included the endpoints (424 and 492). If we try to mimic 
more closely the normal estimate by subtracting one-half the probability at the 
endpoints, we obtain 0.025, essentially the normal number we found earlier. The 
dbinom function Ref. 3 in the R package has been used for this computation. 
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3.3 Shape Comparison Statistics 



There are many different possible statistics for comparing tlie sliapes of tlie 
histograms. We investigate several choices. Table I summarizes the result of each of 
these tests applied to the example in Fig. 1. We list the statistics here, and discuss 
performance in the following sections. 



Even though we don't expect it to follow a distribution, we may evaluate the 
test statistic: 



li Ui — Vi — 0, the contribution to the sum from that bin is zero. We have already 
discussed application of this statistic to the example of Fig. 1. 



Another test statistic we could try may be motivated from a geometric perspec- 
tive. We consider the bin contents of a histogram to define a vector in a A;-dimensional 
space. If two such vectors arc drawn from the same distribution (the null hypothe- 
sis), then they will tend to point in the same direction (we are not interested in the 
lengths of the vectors here). Thus, if we represent each histogram as a unit vector 
with components: 



This is known as the "Bhattacharyya distance measure". We'll refer to it as the 
"BDM" statistic for short. We assume that neither histogram is empty for this 
statistic. All vectors lie in the positive direction in all coordniates, so there is no 
issue with taking the square root. 



3.3.1 Chi-square test for shape 




3.3.2 Geometric test for shape 



{ui/Nu, Uk/Ny}, and {vi/Ny, Vk/N^}, 



we may form the test statistic: 
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It may be noticed that this statistic is related to the statistic - the ■ 
dot product is close to the cross term in the expression. 

We apply this formalism to the example in Fig. 1. The resulting terms in the 
sum over bins are shown in Fig. 4. The sum over bins gives 0.986 (See Table I for 
a summary). According to our estimated distribution of this statistic under the null 
hypothesis, this gives a P-value of 0.97, similar to the test result. 



o 
o 




10 20 30 40 0.95 0.96 0.97 0.98 0.99 

Bin index BDIVI statistic 

Fig. 4. Left: Bin-by-bin contributions to the geometric ("BDM") test statistic for 
the example of Fig. 1. Right: Estimated distribution of the BDM statistic for the 
null hypothesis in the example of Fig. 1. 

3.3.3 Kolmogorov-Smirnov test 

Another approach to a shape test may be based on the Kolmogorov-Smirnov 
(KS) idea. Recall that the idea of the KS test is to estimate the maximum differ- 
ence between observed and predicted cumulative distribution functions (CDFs) and 
compare with expectations. We may adapt this idea to the present case. It should 
be remarked that if we have the actual data points from which the histograms are 
derived, then we may use the Kolmogorov-Smirnov ("KS") procedure directly on 
those points. This would incorporate additional information and yield a potentially 
more powerful test. However, if the bin widths are small compared with possible 
structure it may be expected to not make much difference. 

13 



We modify the KS statistic to apply to comparison of histograms as follows. 
We assume that neither histogram is empty. Form the "cumulative distribution 
histograms" according to: 



Ur 



Then compute the test statistic: 

^KS 



i 



max \Uci 



Test statistics may also be formed for one-tail tests, but we consider only the two-tail 
test here. 

We apply this formalism to the example in Fig. 1. The bin-by-bin distances are 
shown in Fig. 5. The maximum over bins gives 0.043 (See Table I for a summary). 
According to our estimated distribution of this statistic under the null hypothesis, 
this gives a P- value of 0.61, somewhat smaller than for the test result, but still 
indicating consistency of the two histograms. Note that the KS test tends to empha- 
size the region near the peak of the distribution, that is the region where the largest 
fluctuations are expected in Poisson statistics. 
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Fig. 5. Left: Bin- by-bin distances for the Kolmogorov-Smirnov test statistic for 
the example of Fig. 1. Right: Estimated distribution of the Kolmogorov-Smirnov 
distance for the null hypothesis in the example of Fig. 1. 
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3.3.4 Cramer-von-Mises test 



Somewhat similar to the Kolmogorov-Smirnov test is the Cramer-von-Mises 
(CVM) test. The idea in this test is to add up the squared differences between 
the cumulative distributions being compared. Again, this test is usually thought of 
as a test to compare an observed distribution with a presumed parent continuous 
probability distribution. However, the algorithm can be adapted to the two-sample 
comparison, and to the case of comparing two histograms. 

The test statistic for comparing the two samples xi, X2, ■ ■ ■ , xn and yi, y2, ■ ■ ■ , yM 



where is the empirical cumulative distribution for sampling x. That is, Ex{x) — 
n/N if n of the sampled Xi are less than or equal to x. 

We adapt this for the present application of comparing histograms with bin 
contents ui,U2, ■ ■ ■ ,Uk and vi,V2, ■ ■ ■ ,Vk with identical bin boundaries: Let 2; be a 
point in bin i, and define the empirical cumulative distribution function for histogram 
u as: 



We apply this formalism to the example in Fig. 1, finding TcvM = 0.132. The 
resulting estimated distribution under the null hypothesis is shown in Fig. 6. Ac- 
cording to our estimated distribution of this statistic under the null hypothesis, this 
gives a P- value of 0.45 (See Table I for a summary), somewhat smaller than the 
test result. 



is [4]: 




Then the test statistic is: 
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Fig. 6. Estimated distributions of the test statistic for the null hypothesis in the ex- 
ample of Fig. 1 Left: The Cramer-von-Mises statistic. Right: The Anderson-Darling 
statistic. 

3.3.5 Anderson-Dcirling test for shape 

The Anderson-Darling test is another variant on the theme of non-parametric 
comparison of cumulative distributions. It is similar to the Cramer-von-Mises statis- 
tic, but is designed to be sensitive to the tails of the CDF. The original statistic was, 
once again, designed to compare a dataset drawn from a continuous distribution, 
with CDF Fq{x) under the null hypothesis: 

where Fm{x) is the empirical CDF of dataset xi,...Xm- Scholz and Stephens [5] 
provide a form of this statistic for a /c-sample test on grouped data (e.g., as might 
be used to compare k histograms). Based on the result expressed in their Eq. 6, the 
expression of interest to us for two histograms is: 

^ E,(iV„+k-E,) { le^" + - ^"^^1' 
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where k^m is the first bin where either histogram has non-zero counts, kmax is the 
number of bins counting up the the last bin where either histogram has non-zero 
counts, and 

We apply this formahsm to the example in Fig. 1. The resulting estimated 
distribution under the null hypothesis is shown in Fig. 6. The sum over bins gives 
0.849 (See Table I for a summary). According to our estimated distribution of this 
statistic under the null hypothesis, this gives a P-value of 0.45, somewhat smaller 
than the test result, but similar with the CVM result. 

3.3.6 Likelihood ratio test for shape 

We may base a test whether the histograms are sampled from the same shape 
distribution on the same binomial idea as we used for the normalization test. In this 
case, however, there is a binomial associated with each bin of the histogram. We 
start with the null hypothesis, that the two histograms are sampled from the joint 
distribution: 

k Ui Vi 
t—1 

where Ui — a/ii for i — 1,2, . . . , k. That is, the "shapes" of the two histograms are 
the same, although the total contents may differ. 

With ti = Ui + Vi, and fixing the U at the observed values, we have the multi- 
binomial form: 

Pivlu + v = t) = (^^Y^ ( V'"" 

The null hypothesis is that Ui — ajii for all values of i. We would like to test this, 
but there are now two complications: 
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i=l 

3 

= ^^'t'i, and 

i=l 

3 

= ^ ^ ti — ^uj ~l~ ^vj ■ 



i=l 



1. The value of "a" is not specified; 

2. We still have a multivariate distribution. 

For a, we will substitute an estimate from the data, namely the maximum like- 
lihood estimator: 

. Ny 

Nu 

Note that this estimate is a random variable; its use will reduce the effective number 
of degrees of freedom by one. 

We propose to use a likehhood ratio statistic to reduce the problem to a single 
variable. This will be the likelihood under the null hypothesis (with a given by its 
maximum likelihood estimator) , divided by the maximum of the likelihood under the 
alternative hypothesis. Thus, we form the ratio: 



maxjj^ jC{a\v; u + v — t) 
maxHi = ^i/fJ^i}\v;u + v = t) 

n vi+ay \i+a) 



The maximum likelihood estimator, under i^i, for Oj is just 



= Vi/ui. 



Thus, we rewrite our test statistic according to: 



^ _ / l+Vi/Ui f Ny Ui ^ 



\y U,j 



In practice, we'll work with 



i=l 



^ , , 1 + Vi/ui \ , f NyUi 



Before attempting to apply this, we investigate how to handle zero bin contents. 
It is possible that Ui — Vi — for some bin. In this case, P{vi\ui + Vi — ti) — 1, 
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under both Hq and Hi, and this bin contributes zero to the sum. It is also possible 
that ti ^ 0, but Vi — ^ ox Ui — 0. If vi — 0, then 



P(0|ti) 



Vi + A«i 



ti 



Under i^o, this is 



and under H\ it is 



1 + a 



The maximum likelihood estimator for is = 0. Thus, the likelihood ratio for 
bin i is 

ti 



1 + a 

and this contributes to the sum an amount: 

-21nAj = -2tiln 

If instead Ui = 0, then 



Nu + Ny 



P{ti\ti) 

and the contribution to the sum is 



+ [ii 



-21nAi = -2ijln 



We apply this formalism to the example in Fig. 1. The resulting terms in the 
sum over bins are shown in Fig. 7. The sum over bins gives 25.3 (See Table I 
for a summary). This statistic should asymptotically be distributed according to a 

distribution with the number of degrees of freedom equal to one less than the 
number of bins, or A^dof = 39 in this case. If valid, this gives a P-value of 0.96 
in this example. This may be compared with a probability of 0.96 according to the 
estimated actual distribution. In this example we obtain nearly the same answer as 
the naive application of the chi-square calculation with no bins combined. 
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Fig. 7. Value of — 21nAj or Xi as a function of histogram bin in the comparison of 
the two distributions of Fig. 1. Blue circles are — 21nAi; red squares are xf- 

We may see that this close agreement is a result of nearly bin-by-bin equality of 
the two statistics, see Fig. 7. To investigate when this might hold more generally, we 
compare the values of —2 In and xf as a function of Ui and Vi, Fig. 8. We observe 
that the two statistics agree when Ui = Vi with increasing difference away from that 
point. This observation is readily verified analytically. This agreement holds even 
for low statistics. However, we shouldn't conclude that the chi-square approximation 
may be used for low statistics - fluctuations away from equal numbers lead to quite 
different results when we get into the tails at low statistics. Our example doesn't 
really sample these tails. 

The precise value of the probability should not be taken too seriously, except 
to conclude that the two distributions are consistent according to these tests. For 
example, when we combine bins to improve expected x^ behavior, we see fairly large 
fluctuations in the probability estimate just due to the re-binning (Fig. 3). 
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Fig. 8. Value of — 21nAi or Xi as a function of Ui and Vi bin contents. This plot 
assumes Ny^ = Ny. The i subscript is dropped, with the understanding that this 
comparison is for a single bin. 

3.3.7 Likelihood value test for shape 

An often-used but controversial goodness-of-fit statistic is the value of the like- 
lihood at its maximum value under the null hypothesis. It can be demonstrated 
that this statistic carries little or no information in some situations. However, in the 
limit of large statistics it is essentially the chi-square statistic, so there are known 
situations were it is a plausible statistic to use. We thus look at it here. 

Using the results in the previous section, the test statistic is: 
T=-\nC 

If either A^^^ = or Ny 

We apply this formalism to the example in Fig. 1. The resulting estimated 
distribution under the null hypothesis is shown in Fig. 9. The sum over bins gives 
90 (See Table I for a summary). According to our estimated distribution of this 
statistic under the null hypothesis, this gives a P-value of 0.29, similar to the test 
result. The fact that it is similar may be expected from the fact that our example is 
reasonably well-approximated by the large statistics limit. 
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Fig. 9. Estimated distribution of the InjC test statistic for the null hypothesis in the 
example of Fig. 1. 

There are many other possible tests that could be considered, for example, 
schemes that "partition" the to select sensitivity to different characteristics [6] . 

3.4 Distributions Under the Null Hypothesis 

For the situation where the asymptotic distribution may not be good enough, 
we would like to know the probability distribution of our test statistic under the null 
hypothesis. However, we encounter a difficulty: our null hypothesis is not completely 
specified! The problem is that the distribution depends on the values of Ui — a/ii. 
Our null hypothesis only says Vi — a/ii, but says nothing about what /ij might be. 
Note that it also doesn't specify a, but we have already discussed that complication, 
which appears manageable (although in extreme situations one might need to check 
for dependence on a). 

We turn once again to the data to make an estimate for /li, to be used in esti- 
mating the distribution of our test statistics. The straightforward approach is to use 
the maximum likelihood parameter estimators (under Hq): 

fJ-i = ——{ui + Vi), 
I + a 

n = -T-—:[ui^Vi), 
1 + a 
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where a — Ny/Nu- The data is then repeatedly simulated using these values for 
the parameters of the sampling distribution. For each simulation, a value of the 
test statistic is obtained. The distribution so obtained is then an estimate of the 
distribution of the test statistic under the null hypothesis, and P-values may be 
computed from this. Variations in the estimates for /tj and a may be used to check 
robustness of the probability estimates obtained in this way. 

We have just described the approach that was used to compute the estimated 
probabilities for the example of Fig. 1. The bin contents in this case are reasonably 
large, and this approach works well enough for this case. 

Unfortunately, this approach does very poorly in the low-statistics realm. We 
consider a simple test case: Suppose our data is sampled from a flat distribution 
with a mean of 1 count in each of 100 bins. We test how well our estimated null 
hypothesis works for any given test statistic, T, as follows: 

1. Generate a pair of histograms according to the distribution just described. 

(a) Compute T for this pair of histograms. 

(b) Given the pair of histograms, compute the estimated null hypothesis accord- 
ing to the spcified prescription above. 

(c) Generate many pairs of histograms according to the estimated null hypothesis 
in order to obtain an estimated distribution for T. 

(d) Using the estimated distribution for T, determine the estimated P- value for 
the value of T found in step la. 

2. Repeat step 1 many times and make a histogram of the estimated P-values. 
Note that this histogram should be uniform if the estimated P-values are good 
estimates. 
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The distributions of the estimated probabihties for the seven test statistics under 
the null hypothesis are shown in the second column of Fig. 10. If the null hypothesis 
were to be rejected at the estimated 0.01 probability, this algortihm would actually 
reject Ho 19% of the time for the statistic, 16% of the time for the BDM statistic, 
24% of the time for the In A statistic, and 29% of the time for the C statistics, all 
unacceptably larger than the desired 1%. The KS, CVM, and AD statistics are all 
consistent with the desired 1%. For comparison, the first column of Fig. 10 shows 
the distribution for a "large statistics" case, where sampling is from histograms with 
a mean of 100 counts in each bin. We find that all test statistics display the desired 
flat distribution in this case. Table II summarizes these results. 

Table II. Probability that the null hypothesis will be rejected with a cut at 1% on 
the estimated distribution (see text) . Hq is estimated with the bin-by-bin algorithm 
in the first two columns, by the uniform histogram algorithm in the third column, 
and with a Gaussian kernel estimation in the fourth column. 

Test statistic Probability (%) Probabihty (%) Probabihty (%) Probability (%) 
Bin mean = 100 1 1 (uniform) 1 (kernel) 



Hq estimate bin-by-bin bin-by-bin uniform kernel 





0.97 


± 


0.24 


18.5 


± 


1.0 


1.2 


± 


0.3 


1.33 ±0.28 


BDM 


0.91 


± 


0.23 


16.4 


± 


0.9 


0.30 


± 


0.14 


0.79 ±0.22 


KS 


1.12 


± 


0.26 


0.97 


± 


0.24 


1.0 


± 


0.2 


1.21 ±0.27 


CVM 


1.09 


± 


0.26 


0.85 


± 


0.23 


0.8 


± 


0.2 


1.27 ±0.28 


AD 


1.15 


± 


0.26 


0.85 


± 


0.23 


1.0 


± 


0.2 


1.39 ±0.29 


In A 


0.97 


± 


0.24 


24.2 


± 


1.1 


1.5 


± 


0.3 


2.0 ±0.34 


ln£ 


0.97 


± 


0.24 


28.5 


± 


1.1 


0.0 


± 


0.0 


0.061 ±0.061 



24 



Fig. 10. (Figure on previous page) Distribution of the estimated probability that 
the test statistic is worse than that observed, for seven different test statistics. The 
data are generated according to the null hypothesis, consisting of 100 bin histograms 
with a mean of 100 counts (left column) or one count (other columns). The first 
and second columns are for an estimated HO computed as the weighted bin-by-bin 
average. The third column is for an estimated HO where each bin is the average of 
the total contents of both histograms, divided by the number of bins. The rightmost 
column is for an estimated HO estimated with a Gaussian kernel estimator using the 
contents of both histograms. The is computed without combining bins. 

It may be noted that the issue really is one appearing at low statistics. We 
can give some intuition for the observed effect. Consider the likely scenario at low 
statistics that some bins will have zero counts in both histograms. In this case our 
algorithm for the estimated null hypothesis yields a zero mean for these bins. The 
simulation used to determine the probability distribution for the test statistic will 
always have zero counts in these bins, that is, there will always be agreement between 
the two histograms in these bins. Thus, the simulation will find that values of the 
test statistic are more probable than it should. 

If we tried the same study with, say, a mean of 100 counts per bin, we would 
find that the probability estimates are vahd, at least this far into the tails. The left 
column of Fig. 10 shows that more sensible behavior is achieved with larger statistics. 
The x^. In A, and ln>C statistics perform essentially identically at high statistics, as 
expected, since in the normal approximation they are equivalent. 

The AD, CVM, and KS tests are more robust under our estimates of i^o than the 
others, as they tend to emphasize the largest differences and are not so sensitive to 
bins that always agree. For these statistics, we see that our procedure for estimating 
i^o does well even for low statistics, although we caution again that we are not 
examining the far tails of the distribution. 

There are various possible approaches to salvaging the situation in the low statis- 
tics regime. Perhaps the simplest is to rely on the typically valid assumption that 
the underlying distribution is "smooth". Then instead of having an unknown 
parameter for each bin, we only need to estimate a few parameters to describe the 
smooth distribution, and effectively more statistics are available. 
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For example, we may repeat the algorithm for our example of a mean of one 
count per bin, but now assuming a smooth background represented by a uniform 
distribution. This is cheating a bit, since we perhaps aren't supposed to know that 
this is really what we are sampling from, but we'll pretend that we looked at the 
data and decided that this was plausible. As usual, we would in practice want to try 
other possibilities to evaluate systematic effects. 

Thus, we estimate: 

= Nu/k, i = 1,2, 

= Ny/k, i ^ 1,2, . . . ,k. 

The resulting distributions for the estimated probabilities are shown in the third 
column of Fig. 10. These distributions are much more reasonable, at least at the 
level of a per cent (1650 sample experiments are generated in each case, and the 
estimated P value is estimated for each experiment with 1650 evaluations of the null 

hypothesis for that experiment). 

It should be remarked that the ln>C and, perhaps, to a much lesser extent the 
BDM statistic, do not give the desired 1% result, but now err on the "conservative" 
side. It may be possible to mitigate this with a different algorithm, but this has not 
been investigated. We may expect the power of these statistics to suffer under the 
approach taken here. 

Since we aren't supposed to know that our null distribution is uniform, we also 
try another approach to get a feeling for whether we can really do a legitimate 
analysis. Thus, we try a kernel estimator for the null distribution, using the sum of 
the observed histograms as input. In this case, we have chosen a Gaussian kernel, 
with a standard deviation of 2. The "density" package in R [1] is used for this. An 
example of such a kernel estimated distribution is shown in Fig. 11. The resulting 
estimated probability distributions of our test statistics are shown in the rightmost 
column of Fig. 10. In general, this works pretty well. The bandwidth was chosen 
here to be rather small; a larger bandwidth would presumably improve the results. 
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Fig. 11. Sample Gaussian kernel density estimate of the null hypothesis (for sampling 
from a true null). 



3.5 Comparison of Power of Tests 

The power depends on what the alternative hypothesis is. Here, we mostly inves- 
tigate adding a Gaussian component on top of a uniform background distribution. 
This choice is motivated by the scenario where one distribution appears to show 
some peaking structure, while the other does not. We also look briefly at a different 
extreme, that of a rapidly varying alternative. 

The data for this study are generated as follows: The background (null distri- 
bution) has a mean of one event per histogram bin. The Gaussian has a mean of 50 
and a standard deviation of 5, in units of bin number. We vary the amplitude of the 
Gaussian and count how often the null hypothesis is rejected at the 1% confidence 
level. The amplitude is measured in percent, for example a 25% Gaussian has a total 
amplitude corresponding to an average of 25% of the total counts in the histogram, 
including the (small) tails extending beyond the histogram boundaries. The Gaus- 
sian counts are added to the counts from the null distribution. An example is shown 
in Fig. 12. 
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Fig. 12. Left: The mean bin contents for a 25% Gaussian on a flat background of 
one count/bin (note the suppressed zero). Right: Example sampling from the 25% 
Gaussian (filled blue dots) and from the uniform background (open red squares). 



The distribution of estimated probability, under HO, that the test statistic is 
worse than that observed (i.e., the distribution of P- values) is shown in Fig. 13 for 
seven different test statistics. Three different magnitudes of the Gaussian amplitude 
are displayed. The power of the tests to reject the null hypothesis at the 99% confi- 
dence level is summarized in Table III and in Fig. 14 for several different alternative 
hypothesis amplitudes. 



Table III. Estimates of power for seven different test statistics, as a function of Hi. 
The comparison histogram {Hq) is generated with all = 100 bins Poisson of mean 
1. The selection is at the 99% confidence level, that is, the null hypothesis is accepted 
with (an estimated) 99% probability if it is true. 





HO 


12.5 


25 


37.5 50 -25 


Statistic 


% 


% 


% 


% % % 




1.2 ±0.3 


1.3 ±0.3 


4.3 ±0.5 


12.2 ±0.8 34.2 ± 1.2 1.6 ± 0.3 


BDM 


0.30 ±0.14 0.5 ±0.2 


2.3 ±0.4 


10.7 ±0.8 40.5 ± 1.2 0.9 ±0.2 


KS 


1.0 ±0.2 


3.6 ±0.5 13.5 ±0.8 


48.3 ± 1.2 91.9 ±0.7 7.2 ± 0.6 


CVM 


0.8 ±0.2 


1.7±0.3 


4.8 ±0.5 


35.2 ± 1.2 90.9 ±0.7 2.7 ±0.4 


AD 


1.0 ±0.2 


1.8 ±0.3 


6.5 ±0.6 


42.1 ± 1.2 94.7 ±0.6 2.8 ± 0.4 


In A 


1.5 ±0.3 


1.9 ±0.3 


6.4 ±0.6 


22.9 ± 1.0 67.1 ±1.2 2.4 ±0.4 


\nC 


0.0 ±0.0 


0.1 ±0.1 


0.8 ±0.2 


6.5 ±0.6 34.8 ±1.2 0.0 ± 0.0 
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Fig. 13. (Figure on previous page) Distribution of estimated probability, under HO, 
that the test statistic is worse than that observed, for seven different test statistics. 
The data are generated according to a uniform distribution, consisting of 100 bin 
histograms with a mean of 1 count, for one histogram, and for the other histogram 
with a uniform distribution plus a Gaussian of strength 12.5% (left column), 25% 
(middle column), and 50% (right column). The is computed without combining 
bins. 




Gaussian Amplitude (%) 

Fig. 14. Summary of power of seven different test statistics, for the alternative 
hypothesis with a Gaussian bump. Left: linear vertical scale; Right: logarithmic 
vertical scale. [Best viewed in color. At an amplitude of 35%, the ordering, from top 
to bottom, of the curves is: KS, AD, CVM, In A, x'^, BDM, ln£.] 



In Table IV we take a look at the performance of our seven statistics for his- 
tograms with large bin contents. It is interesting that in this large-statistics case, for 
the and similar tests, the power to reject a dip is greater than the power to reject 
a bump of the same area. This is presumably because the "error estimates" for the 
are based on the square root of the observed counts, and hence give smaller errors 
for smaller bin contents. We also observe that the comparative strength of the KS, 
CVM, and AD tests versus the x^, BDM, In A, and ln£ tests in the small statistics 
situation is largely reversed in the large statistics case. 
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Table IV. Estimates of power for seven different test statistics, as a function of Hi. 
The comparison histogram (Hq) is generated with all /c = 100 bins Poisson of mean 
100. The selection is at the 99% confidence level. 





HO 5 


-5 


btatistic 


% % 


% 


X 


0.91 ± 0.23 79.9 ± 1.0 92.1 ± 0.7 


EDM 


0.97 ±0.24 80.1 ± 1.0 


92.2 ±0.7 


KS 


1.03 ±0.25 77.3 ± 1.0 


77.6 ± 1.0 


CVM 


0.91 ±0.23 69.0 ± 1.1 


62.4 ± 1.2 


AD 


0.91 ±0.23 67.5 ± 1.2 


57.8 ± 1.2 


In A 


0.91 ±0.23 79.9 ± 1.0 


92.1 ±0.7 


ln£ 


0.97 ±0.24 79.9 ± 1.0 


91.9 ±0.7 



To get an idea of what happens for a radically different alternative to the null 
distribution, we consider sensitivity to sampling from the "sawtooth" distribution as 
shown in figure 15. This is to be compared once again to samplings from the uniform 
histogram. The results are tabulated in Table V. The "percentage" sawtooth here 
refers to the fraction of the null hypothesis mean. That is, a 100% sawtooth on a 1 
count/bin background oscillates between a mean of counts/bin and 2 counts/bin. 
The period of the sawtooth is always two bins. 
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Fig. 15. Left: The mean bin contents for a 50% sawtooth on a flat background 
of one count/bin (blue), compared with the fiat background means (red). Right: 
Example sampling from the 50% sawtooth (filled blue dots) and from the uniform 
background (open red squares). 
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In this example, the and hkehhood ratio tests are the clear winners, with 
BDM next. The KS, CVM, and AD tests reject the null hypothesis with the same 
probability as for sampling from a true null distribution. This very poor performance 

for these tests is readily understood, as these tests are all based on the cumulative 
distributions, which average out local oscillations. 

Table V. Estimates of power for seven different test statistics, for a "sawtooth" 
alternative distribution. 





50 


100 


Statistic 


% 


% 




3.7 ±0.5 


47.8 ± 1.2 


BDM 


1.9 ±0.3 


33.6 ± 1.2 


KS 


0.85 ±0.23 


1.0 ±0.2 


CVM 


0.91 ±0.23 


1.0 ±0.2 


AD 


0.91 ±0.23 


1.2 ±0.3 


In A 


4.5 ±0.5 


49.6 ± 1.2 


ln£ 


0.30 ±0.14 10.0 ±0.7 



4. Conclusions 

These studies have demonstrated some important lessons in "goodness-of-fit" 
testing: 

1. There is no single "best" test for all apphcations. Statements such as "test X 
is better than test Y" are empty without giving more context. For example, 
the Anderson-Darling test is often very powerful in testing normality of data 
against alternatives with non- normal tails (such as the Cauchy distribution) [7] . 
However, wc have seen that it is not always especially powerful in other situations. 
The more wc know about what we wish to test for, the more reliably wc can 
choose a powerful test. Each of the tests investigated here may be reasonable 
to use, depending on the circumstance. Even the controversial C test works as 
well as the others sometimes. However, there is no known situation where it 
actually performs better than all of the others, and indeed the situations where 
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it is observed to perform as well are here limited to those where it is equivalent 
to another test. 

2. Computing probabilities via simulations is a very useful technique. However, it 
must be done with care. The issue of tests with an incompletely specified null 
hypothesis is particularly insidious. Simply generating a distribution according 
to some assumed null distribution can lead to badly wrong results. Where this 
could occur, it is important to verify the validity of the procedure. Note that we 
have only looked into the tails to the 1% level. The validity must be checked to 
whatever level of probability is needed for the results. Thus, we cannot blindly 
assume the results quoted here at the 1% level will still be true at, say, the 0.1% 
level. 

We have concentrated in this paper on the specific question of comparing two 
histograms. However, the general considerations apply more generally, to testing 
whether two datasets are consistent with being drawn from the same distribution, 
and to testing whether a dataset is consistent with a predicted distribution. The 
KS, CVM, AD, ln£, and L tests may all be constructed for these other situations 
(as well as the and EDM, if we bin the data). 
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